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A  BSTBACT 


This  thesis  uses  renewal  theory  to  investigate  the 
lanchester-type  combat  attrition  process.  The  attrition 
process  is  analysed  in  detail  and  modelled  as  a  so-called 
renewal  process  in  which  times  between  casualties  are 
considered  to  be  independently  and  identically  distributed 
random  variables.  Other  random  variables  that  can  be 

considered  in  the  renewal  process  are  examined,  and  the 
distributions  of  these  random  variables  are  determined  in 
order  to  study  the  behaviour  of  attrition  process.  Examples 
with  specific  distribution  functions  are  given  for  better 
understanding.  Computer  simulation  is  generated  and  compared 
with  the  attrition  model  developed.  The  total  casualty 
occurence  by  total  force  is  also  discussed,  through  pooling 
of  the  single  renewal  processes.  The  total  casualty  occu¬ 
rence  is  shown  to  be  a  Poisson  process  and  times  between 
casualties  to  be  approximately  exponentially  distributed  for 
large  numbers  of  combatants. 


4 


TABLE  OF  CONTESTS 


I.  INTRCDUCT  ION . 9 

II.  FUNDAMENTAL  L INCHES TER- TYPE  ATTRITION  PARADIGM  .  .  12 

A.  THE  BASIC  PARADIGM . 12 

B.  DETERMINATION  OF  EXPECTED  TIME  TO  KILL  ....  14 

C.  TARGET  ACQUISITION  PROCESS  .  15 

1.  Serial  Acquisition . 16 

2.  Parallel  Acquisition . ..17 

III.  ANALYSIS  OF  THE  COMBAT  ATTRITION  PROCESS  .  20 

A.  INTRODUCTION . 20 

B.  REVIEW  OF  RENEWAL  PROCESS . 22 

C.  DISTRIBUTION  OF  TOTAL  NUMBER  OF  CASUALTIES  .  .  25 

D.  EXPECTED  NUMBER  OF  CASUALTIES  .  26 

1.  Renewal  Function  . . 26 

2.  Renewal  Equation  . . ...30 

3.  Renewal  Density  . . 31 

E.  LIMIT  BEHAVIOUR  OF  CASUALTY  OCCURENCE  ....  32 

F.  AGS  AND  EXCESS  LIFE  DISTRIBUTION  OF 

CASUALTY  OCCURENCE  .  34 

IV.  THE  SUPERPOSITION  OF  CASUALTY  OCCURENCES  .  37 

A.  INTRODUCTION . 37 

B.  SOME  GENERAL  PROPERTIES . 38 

1.  Normal  Approximation  of  Pooled  Output  .  .  38 

2.  Poisson  Process  of  Pooled  Output  .  ....  41 

V.  COMPARISON  OF  RENEWAL  PROCESS  WITH  COMPUTER 

SIMULATION . 45 


\ 

1 


5 


Gyv*  •  ■ 


VI.  CONCLUSIONS . 43 

A.  SUM  M  A  El . 43 

3.  SUGGEST IO MS  FOR  FURIHEfl  STDDX  .  50 

APPENDIX  A:  FORTRAN  PROGRA MMING  FOR  SIMULATING  THE 

KILL  TIME  (N  =  3} . 51 

APPENDIX  B:  SIMULATED  DATA  FOR  KILL  TIME  (N  =  3,  M  » 

6) . 52 

APPENDIX  C:  DISCRETE  APPROXIMATION  FO  RENEWAL 

FUNCTION . 53 

LIST  OF  REFERENCES . 54 

INITIAL  LISTRIBUTION  LIST  .  55 


6 


LIST  OF  TABLES 


I.  Variables  for  Expected  Time  to  Kill . 15 

II.  Attrition  Process  Using  Serial  Acquisition 

Mode . 17 

III.  Attrition  Process  Using  Parallel  Acquisition 

Mode . 18 

IV.  Random  Variables  in  Renewal  Process  ...  ....  24 


LIST  OF  FIGURES 


3.1  Renewal  Process  . . 25 

3.2  Age  and  Excess  Life  of  Casualty  Occurence  ...  35 

4.1  Superposition  of  Casualty  Occurences . 37 

4.2  Superposition  of  Excess  Life . 42 

5.1  Expected  Kills  From  Theoretical  And 

Simulated  . . 43 


8 


I.  INTRODUCTION 


Today#  deterministic-differential-equation  models (  that 
are  commonly  called  Lanchester-type  combat  models)  are 
widely  used  in  military  operations  research  for  defense 
planning  purposes#  especially  as  concerns  land  combat.  In 
combat  analysis#  such  questions  frequently  arise  as  'Who 
will  win  the  battle  ?#  How  long  will  the  battle  last  ?,  How 
many  survivors  will  the  winner  have  ?,  and  How  do  the  force 
levels  change  over  time  ?'#etc.  There  have  been  a  lot  of 
studies  to  estimate  attrition  in  battle  to  answer  such  ques¬ 
tions  as  these. 

In  1914#F.  W.  Lanchester  [Sef.  1]  first  proposed  a  mathe¬ 
matical  model  of  combat  in  order  to  justify  the  principle  of 
concentraticn  of  forces  under  "modern"  conditions.  He 
believed  that  in  "ancient"  times  war  was  a  series  of 
one-to-one  duels  between  men  engaged  in  hand-to-hand  combat 
using  sword,  axe,  lance,  etc.  So  tnere  would  be  no  advan¬ 
tage  gained  from  concentration  of  forces.  Lanchester 
reasoned  that  modern  technology  had  changed  the  nature  of 
this  "ancient"  warfare  and  now  there  was  a  decided  advantage 
to  "concentration  of  forces". i  He  hypothized  that  each 
side's  casualty  rate  was  proportional  to  the  number  of 
opposing  firers.  Implicit  (  but  not  yet  explicit)  in 
Lanchester 's  original  work  is  the  concept  of  an  "attrition- 
rate  coefficient"  as  the  rate  at  which  a  single  firer  kills 
a  particular  enemy  target  type.  Furthermore,  Lanchester 


1  Helmbold  notes  that  the  analysis  of  the  data  from  92 
land  combat  battles  suggests  that  victory  in  cattle  is 
primarily  determined  by  the  factors  other  than  numerical 
superiority  and  challenges  the  ability  of  any  model  of 
combat  which  concentrates  almost  exclusively  on  numerical 
force  size  to  yield  a  practically  useful  predictor  of 
victory  in  battle. 
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tnought  that  the  number  of  survivors  were  the  critical 
factor  deciding  the  outcome  of  the  battle  in  "modern" 
warfare. 

Since  Lanchester's  original  efforts  to  describe  the 
dynamics  of  combat  mathematically,  many  analysts  -  in  an 
attempt  to  add  reality  to  the  combat  description  -  have 
extended  the  theory  to  include  additional  factors.  Many  of 
these  extentions  were  described  by  Bolansky  [Ref.  2].  In 
his  paper,  Dolansky  pointed  out  that  the  Lanchester 
attrition-rate  coefficients  were  aard  to  determine  for 
particular  weapons,  and,  accordingly,  coefficients  of  known 
magnitude  had  been  assumed  in  most  models.  Thus,  Dolansky 
concluded  that  the  use  of  Lanchester-type  models  for 
prediction  of  battle  results  had  been  hampered  by  this 
inability  to  predict  numerical  values  for  attrition-rate 
coefficients.  However,  even  if  we  may  be  in  doubt  as  to  the 
proper  form  of  attrition  relations,  we  may  at  least  be  able 
to  make  pragmatic  assertions  such  as,  *  If  the  Lanchester 
attrition  is  used  for  predictions,  then  by  measuring  certain 
factors  and  performing  certain  operations  a  prediction  of 
such-and-such  an  aspect  of  combat  can  be  obtained. '  So, 
although  there  has  been  a  continuing  discussion  among  nd  li- 
tary  operations  analysts  about  the  merits  of  Lanchester-type 
attrition-rate,  it  is  still  generally  accepted  as  the  basic 
founding  for  all  the  combat  models  studies  and  attrition- 
rate  coefficient  is  still  a  core  of  these  studies. 

Bonder/Barf oot  [Ref.  3]  proposed  thatian  attrition-rate 
coefficient  be  defined  as  the  reciprocal  of  the  expected 
time  for  an  individual  firer  to  kill  an  enemy  target.  A 
more  thorough  dicussion,  however,  of  the  justification  for 
Bonder/Barefoot  is  to  be  found  in  Taylor  [Ref.  4],  Barf  oot 
[Ref.  5]  suggested  that  the  attrition  process  be  considered 
as  a  renewal  process.  Thus,  for  a  singxe  Y  firer,  we  will 
assume  that  each  individual  Y  kills  X  targets  according  to 


an  attrition  process  an  which  tames  to  kills  are  indepen¬ 
dently  and  identically  distributed  (i.i.d.)  random  variables 
(r.v).  He  will  also  assume  that  each  firer  acts  indepen¬ 
dently  ox  any  other  firer,  and  similarly  for  targets.  Suca 
an  attrition  process  turns  out  to  be  a  renewal  process  in 
the  theory  of  stochastic  processes. 

Once  we  introduce  the  theory  of  renewal  process,  we  car. 

consider  many  other  random  variables  such  as  the  time  up  to 

n-th  kill,  S  ;  the  number  of  kills  in  time  t,  !J  (t)  ;  the 
n 

expected  number  of  kills  in  time  t,  :i(t),  etc.  Thus  this 
thesis  starts  with  consideration  of  3onder/Taylor 's  concept 
of  attrition-rate  coefficient,  and  then  studies  the  prob¬ 
ability  distribution  of  the  time  between  casualties 
according  to  the  theory  of  renewal  process. 

In  this  thesis,  we  will  mainly  focus  on  the  fundamental 
Lanchester-type  attrition  paradigm  waich  deals  with  homoge¬ 
neous  forces,  because  understanding  this  basic  paradigm  is 
essential  fcr  extenticns  to  further  complex  medals. 

This  thesis  is  organized  in  the  following  fashion. 
First,  we  review  a  few  relevant  facts  about  the  Lanchester 
attrition  models.  Then  we  review  oasic  ideas  of  renewal 
process  and  develop  the  random  variables  whose  study  is  of 
our  concern,  followed  by  the  study  of  these  random  vari¬ 
ables.  Then  we  will  describe  the  combat  attrition  process  by 
comparing  the  values  of  different  random  variables  obtained 
through  renewal  process.  Fanally,  the  conclusion  mainly 
contributes  to  the  summary  of  the  work  and  suggests  further 
studies. 


II.  FONDA MENIAL  LAMCHB S TER- TYPE  ATTRITION  PARADIGM 

A.  THE  EASIC  PARADIGM 

Let  us  consider  combat  between  two  homogeneous  forces:  a 
homogeneous  X  force  opposed  by  a  homogeneous  Y  force.  We 
will  focus  on  the  f orce-on -force  attrition  process  in  the 
combat  between  these  two  homogeneous  forces.  The  basic 
Lanchester-type  paradigm  for  "modern"  warfare  assumes  that 
the  casualty  rate  of  such  a  homogeneous  force  is  directly 
proportional  to  the  number  of  enemy  forces,  e.g.  the  X  force 
casualty  rate  is  given  by 

dx 

--  =  -  ay  with  x(3)  =  x  , 

dt  0 

and  the  Y  force  casualty  rate  is  given  by 

dy 

--  *  -  bx  with  y(3)  »  y  , 

dt  0 

where  *a'  denotes  the  rate  at  which  a  single  Y  hirer  kills  X 

targets  and  is  called  a  Lan Chester  attrition- rate  coef f i- 

cient.  Here  x  and  y  denote  the  number  of  X  and  Y 
-  0  0 

comcatants  (respectively)  at  the  beginning  of  the  battle. 
Many  factors  may  affect  the  value  of  'a’  and  'b',  but  for 
tae  time  being  it  is  not  essential  that  we  be  explicit  about 
functional  dependence  of  a  and  b.  Let  us  next  address  the 
problem  of  computing  a  reliable  numerical  values  for  such 
attrition-rate  coefficients. 

Two  approaches  have  been  developed  in  the  United  States 
to  determine  the  value  of  these  coefficients  as  follows: 

(1)  A  statistical  estimate  based  on  "combat"  data 
generated  by  a  detailed  Monte  Carlo  combat  simulation. 
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(2)  An  analytical  submodel  or  attrition  process  for  tie 
particular  combination  of  firer  and  target  type. 

In  this  thesis  we  will  consider  only  tne  second  approach. 
S. Bonder  [Ref.  6]  has  called  this  approach  the  use  of  a 
free-standing  or  independent  analytical  mod el.  The  basic 
idea  of  this  approach  is  to  develop  an  analytical  expression 
for  each  required  kill  rate  by  considering  the  single  firer 
engaging  a  'passive'  target  (i.e.  one  that  doesn't  fire 
back)  and  then  tie  ali  the  attrition  rates  together  in 
force-on-force  combat  with  a  Lanche ster-type  model. 

Barfoot  [Bef.  5]  in  his  study  has  suggested  taking  a 
Lanchester  attrition-rate  coefficient  as  the  reciprocal  of 
the  expected  time  for  an  i ndi vidual  firer  to  kill  an  enemy 
'target. '2  flithin  ccntext  of  the  above  homogeneous  force 
Lanchester-type  combat  model, this  means 


a  = 


1 

2[T 


XI 


where  T  denotes  a  r.v.  representing  the  time  for  a  single 
XY 

I  firer  to  kill  an  X  target  and  Z[T]  denotes  expected  time 
to  kill.  Taylor  also  supported  Barfoot’s  idea  in  his  recent 
study  [Ref.  7]  by  comparing  times  between  kills  that  were 
both  exponentially  distributed  and  non-exponentially 
distributed,  and  using  previous  equation  to  determine  the 
expected  time  to  kill  a  target. 


^However,  the  justification  of  this  approach  is  npt 
accepted  and  xs  apparently  somewhat  controversial.  But  m 
case  of  exponentially  distributed  times  between  kills,  a 
deterministic  Lanchester-type  model  may  indeed  be  considered 
to  yield  the  mean  course  of  combat.  Anyway,  this  approach  is 
accepted  as  general  guide  for  further  study  of  attrition- 
rate  coefficient. 


B.  OETEB 91  NATION  OF  EXPECTED  TIME  TO  KILL 

Bonder  and  Farrel  ££ef.  8]  nave  developed  general  meth¬ 
odology  for  determining  the  attrition-rate  coefficient  by 
determining  the  expected  time  to  itill  a.  target,  2[T],  for  a 
vide  spectrum  of  weapon  system  types.  Morever,  research 
since  the  aid-1960's  has  led  to  the  development  of  several 
other  methods  for  computing  the  expected  time  to  kill  a 
target.  For  present  purposes,  two  methods  will  he  of  our 
interest, 

(1)  Method  based  on  sum  of  component  event  times, 
and 

(2)  Method  based  on  first-passage  time  in  Semi-Markov 
process 

If  we  assume  the  special  case  of  tactical  interest:  namely, 
tne  case  of  Markov-detenden  t  fire.  Bonier  [Ref.  3]  has  shown 
that  E[  T  ]  turns  out  to  ne  as  follows: 


(t  *t  )  (t  *t  ) 

n  f  m  f  r  1-?  (fl|  a)  ] 

E£T  ]=t  +  t  -t  + - - {- - -  ♦  P  ( H  |  H)  -  ?  } 

a  1  a  P  <K  1 H)  ?(9M)  ?(K|H)  1 


where  all  symbols  are  defined  in  Table  I  This  expession  for 
E[T]  holds  for  the  fcilowing  assumptions; 

(1)  Markov-dependent  fire  with  parameter  ?  ,  P  (H|  H)  , 
and  ?  (H|M) 

(2)  Geometric  distribution  for  tne  number  of  hits 
required  for  a  kill  with  parameter  P(K|H) 

In  assessing  the  time  factor  to  determine  the 
attrition-rate  coefficient,  it  turned  out  that  target  acqui¬ 
sition  process  and  weapon  system  capanility  were  most  impor¬ 
tant  factor  from  the  formula.  However,  when  we  consider  the 
time  factor  of  the  attrition  process,  target  acquisition 
process  would  be  the  most  important  factor  to  study  because 
it  directly  affects  the  required  length  of  time  to  kill  a 
target.  So  we  will  elaborate  on  target  acquisition  process 
for  further  use  in  attrition  analysis. 
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TABLE  I 

Variables  for  Expected  Tine  to  Kill 


”1 


time  to  acquire  a  target. 


t 


a 

time  to  fire  first  round  after  target  acguired, 
time  to  fire  a  round  following  a  hit,t 

h 

time  to  fire  a  round  following  a  miss,t 


t 


1 


m 


time  of  flight 
probability  of 
probability  of 
probability  of 
probability  of 


of  the  projectile,  t _ 
a  hit  on  first  round, 

a  hit  on  a  round  following  a  hit,  P(H|H) 
a  hit  on  a  round  following  a  miss,  ?  (K|  M) 
destroying  a  target  given  hit,  P(KjH) 


j 


C.  TAB  GET  ACQUISITION  PROCESS 

Target  acguisition  is  the  initial  step  in  engaging  a 
target  and  a  most  deciding  factor  that  affects  the  time  to 
kill.  An  important  distinction  made  in  VEI’s  Lanchester-type 
combat  models  is  whether  the  target  acguisition  process  of  a 
single  "typical"  firer  type  is  a  serial  process  or  a 
parallel  process  [Bef.  4].  The  two  modes  for  the  target 
acguisition  process  considered  by  VRl's  models (  including 
VSCTOR-2  )  are  as  follows: 

(1)  Serial  Acguisition 

(2)  Parallel  Acguisition 
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Here  a  single  firer  using  serial  acquisition  does  not 
acquire  targets  while  engaging  another  target,  say  target  A. 
When  such  a  rirer  ceases  to  engage  a  target  A  (  due  to  hill, 
lost  or  any  other  reason),  then  he  must  acquire  a  new  target 
all  the  way  from  the  beginning.  So  it  is  assumed  that  he 
does  not  remember  any  acquisitions  made  prior  to  engaging 
the  target  A.  On  the  other  hand,  a  firer  using  parallel 
acquisition  searchs  for  targets  continuously,  e7en  while 
engaging  target  A  and  remembers  those  targets  that  have  been 
acquired.  Thus  when  such  a  firer  finishes  the  engagement 
with  target  A,  he  can  immediately  snift  engagement  to  the 
next  target,  provided  that  such  a  target  was  acquired  and 
processed  during  or  before  the  engagement  of  target  A.  For 
both  methods,  the  VRI  models  assume  that  a  firer  never 
engages  a  killed  target  again. 

1.  Serial  Acquisition 

Now  when  we  consider  the  basic  Lanchester-type 
combat  model  in  the  previous  section,  a  firer  in  serial 
process  mode  of  target  acquisition  starts  again  from  the 
beginning  in  searching  a  new  target  whenever  the  previous 
target  has  teen  killed.  Thus  we  assume  that  target  acquisi¬ 
tion  time  is  always  needed  to  acquire  a  different  target  for 
engagement.  Here  it  is  assumed  that  the  total-force  kill 
rate  is  just  the  single-fir er-kill-rate  times  the  number  of 
firers,  e.g. 


dx 


since  each  firer  is  independent  from  any  other  firer.  Also, 
in  the  case  of  a  serial  process  for  each  firer,  the  expected 
time  to  kill  a  target  is  the  same  as  has  been  shown  in  the 
equation  for  E[  T  ]  in  section  B.  A  summary  of  attrition-rate 
coefficient  results  for  serial  acquisition  is  given  in 
Table  II  . 
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TABLE  II 

Attrition  Process  Using  Serial  Acguisition  Mode 


dx 

—  =  -  ay 
dt 


2£I  ] 

XY 


T  =  time  for  a  Y  firer  to  kill  an  X  target 
XY 

(t  *t  )  (t  +t  ) 
h  f  ffl  f  [  1-p  (H  J  h)  ] 

EfT]  *  t  *t  -t  ♦ - + - {- - -  ♦  ?  ( 3 1  H)  -  P  } 

a  1  h  P(K|  H)  P  (HI  M)  ?  (K |  H)  1 


2 «  Parallel  Acquisition 

As  we  aenticned  earlier,  a  firer  in  the  parallel 
mode  of  acquisition  continues  to  acquire  new  targets  even 
while  be  is  engaging  a  particular  target.  Once  this  target 
has  been  killed,  then  he  can  ismediately  shift  fire  to  a  new 
target  which  has  been  acquired  already.  In  this  case, 
total-force  kill  rate  is  given  by  the  product  of  the  kill 
rate  of  single  firer  against  acquired  targets  and  the 
expected  nuafcer  of  firers  who  have  already  acquired  one  or 
sore  targets. 

The  sunaary  results  for  attrition  process  using 

parallel  acguisition  is  given  in  table  III  .  Note  here  that 

f  is  same  as  the  probability  that  a  single  Y  firer  using 
xy 


li 
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parallel  acquisition  node  has  available  one  or  more  acquired 

targets  at  which  to  fire  at  tine  t,  because  if  he  is  firing 

we  know  that  he  already  acquired  the  targets.  Then  it  should 

also  be  noted  that  f  y  represents  the  expected  number  of  Y 

xy 

firers  who  have  already  acquired  one  or  more  X  targets. 


TABLE  III 

Attrition  Process  Using  Parallel  Acquisition  node 


dx 

dt 


f 

xy 


a  *  y 


f  :  Probfsingle  Y  finer  using  parallel  acquisition 
xy  is  firing  at  X  target  at  random  tine  t  } 


't 

f  =  1  -  exp  {  -  x  fi  (s)  ds} 
xy  *  0  x  y 


a* 


1 


E[T«  ] 
X  Y  J 


:Y  firer  kill  rate  against  acquired  target 


I1  =  Tine  for  a  Y  firer  to  kill  an  acquired  X  target 
XY  (conditional  kill  tine) 


E  :  the  rate  at  which  a  Y  firer  acquires  X  targets 
xy  at  tine  t  when  appropriate 


E[T»  ] 


(t  *t  ) 
h  f 

t  -t  «■ - 

1  h  P(KJH) 


♦ 


(t  *t  ) 

a  f  £  1  — P  (H|  H)  1 
- (h— iJ  +  p 

?(H|d)  P(KjH) 


(H  |  H)  } 
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Nov  ve  stay  assume  that  it  wo-uld  take  longer  to  kill 
a  target  for  a  firer  using  serial  acquisition  than  a  firer 
using  parallel  acquisition  because  or  the  time  fraction  of 
acquisition  of  target.  This  distinction  between  these  two 
acquisition  nodes  may  be  further  studied  as  the  case  of  a 
delayed  renewal  process.  Here  we  assume  that  expected  time 
to  kill  the  first  target  would  be  obtained  from  the  serial 
acquisition  but  following  time  to  kill  would  be  from  the 
parallel  acquisition,  provided  that  the  targets  are  acquired 
already.  But  in  this  thesis,  we  will  only  deal  with  the 
firer  who  uses  the  same  acquisition  mode  to  attrit  all 
targets. 

Before  continuing,  it  is  important  that  the  reader 
understand  both  the  basic  Lanchester-type  attrition  model 
and  also  Bonder’s  approach  of  tne  attrition-rate  coefficient 
as  the  reciprocal  of  expected  time  to  kill  a  target.  It  must 
be  recognized  that  this  expected  tine  to  kill  a  target  is 
primarily  dependent  cn  the  acquisition  process. 
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III.  ANALYSIS  QZ  JHE  COMBAT  ATTRITION  PROCESS 


A.  INTRODUCTION 

The  Lanchester  attrition-rate  coefficient  is  the  rate  at 
which  a  single  firer  kills  a  particular  enemy  target  type  in 
Lanchester-type  combat.  Development  of  technically-sound 
and  scienti ficaily- valid  methodology  for  determining  numer¬ 
ical  values  for  Lanchester  attrition-rate  coefficient  is  an 
essential  prerequisite  for  building  a  Lanchester-type  combat 
model.  The  basic  construct  cf  3onder/3arf oot  methodology  was 
to  take  a  Lanchester  attrition-rate  coefficient  as  the 
|  reciprocal  cf  the  expected  time  for  an  individual  firer  to 

kill  an  enemy  target  as  mentioned  in  the  previous  chapter. 
Taylor  provided  justification  [Ref.  4]  for  taking  Lanchester 
attrition-rate  coefficient  as  the  reciprocal  of  the  expected 
time  to  kill.  He  started  from  the  oasic  hypothesis  that 
combat  is  a  complex  random  process,  but  it  contains  enough 
regularity  that  the  appropriate  Lanchester-type  equations 
are  a  gocd  approximation  to  the  mean  course  of  combat.  It  is 
clear  that  the  casualty  rate  is  equal  to  the  reciprocal  of 
the  expected  time  for  a  force  to  inflict  a  casualty,  when 
the  times  between  casualties  are  exponentially  distributed. 
However,  in  case  the  times  between  casualties  are  no  longer 
exponentially  distributed,  Taylor  used  the  suggestion  made 
by  Bonder  ana  Barfoot.  Bonder  [Ref.  3]  and  Barf  oot  [Ref.  5] 
suggested  defining  the  Lancnester  attrition-rate  coefficient 
as  the  expected  rate  which  a  single  firer  kills  enemy 
targets. 

In  the  spirit  of  Bonder  and  Barfoot,  Taylor  provided 
more  rigorous  justification  [Ref.  7]  for  the  Lanchester 
attrition-rate  coefficient  that  does  not  assume  an 
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exponential  distribution  for  tines  between  casualties.  He 
considered  the  case  where  the  initial  force  size  of  X  ana  Y 
is  large  enough  to  insure  a  negligible  probability  that  the 
battle  is  terminated  before  an  attrition  pattern  coula  be 
established.  He  made  no  specific  assumptions  acout  the 
distribution  of  times  between  kills,  but  assumed  that  each 
individual  Y  force  kills  X  targets  according  to  an  attrition 
process  in  which  the  times  between  kills  are  independently 
and  identically  distributed  random  variables.  Thus  in  the 
parlance  of  the  theory  of  stochastic  processes,  he  said  that 
such  an  attrition  process  is  called  a  renewal  process. 
Prior  to  this,  attrition  prediction  nas  been  difficult  due 
to  the  inability  to  predict  casualty  patterns.  By  investi¬ 
gating  the  distribution  of  the  times  between  casualties  via 
renewal  theory,  combat  attrition  analysis  becomes  easier. 
Indeed,  we  know  a  great  deal  about  the  casualty  pattern  once 
we  know  the  distribution  of  the  interarrival  time  of  kills. 
As  we  consider  more  variables  from  the  renewal  process,  we 
find  more  information  about  tne  casualty  pattern. 
Furthermore,  we  can  track  the  development  of  the  combat 
attrition  more  precisely  by  studying  the  probabilistic 
distribution  function  of  those  random  variables  that  can  be 
considered  in  the  renewal  process. 

We  assume  here  that  the  reader  is  aware  of  the  concept 
of  a  counting  process.  Wnen  we  think  about  combat  between  X 
and  Y  forces,  we  can  assume  that  time  between  kills  by  a 
single  firer  has  some  distribution  F.  Now  we  are  interested 
in  observing  the  occurence  of  casualty  and  the  number  of 
casualties,  N(t),  that  have  happened  in  time  interval  (0,t]. 
Furthermore,  when  times  between  kills  are  i.i.d.,  we  call 
this  counting  process  as  renewal  process. 

In  this  chapter,  the  study  of  attrition  analysis  will 
include  the  following  : 


r 

► 


•  The  r.v.,  N  (t)  ,  total  number  of  casualties  by  tire  t, 
it’s  distribution,  and  properties. 

«  The  expected  number  of  casualties  uy  tiae  t. 

•  Theorems  discussing  casualty  occurence. 

•  The  limiting  behavior  of  the  casualty  occurence. 

•  The  age  and  excessive  life  of  the  casualty  occurence. 

•  Examples  with  specific  casualty  distributions . 

•  Application  of  renewal  process  for  casualty  estimation. 

•  Total  force  attrition  by  pooling  the  single  firer 
renewal  processes. 

B.  REVIEW  OF  RENEW  A 1  PROCESS 

If  the  sequence  of  nonnegative  random  variables  {X  ,X  , 

1  2 

...,}  are  i.i.d.,  then  the  counting  process  (H(t),  t>0  }  is 
said  to  be  a  renewal  process  (or  often  called  an  ordinary 
renewal  process) .  Thus,  a  renewal  process  is  a  counting 

process  such  that  the  time  until  the  first  event  has  seme 

distribution  F,  the  time  between  the  first  and  the  second 
event  has,  independently  of  the  tiae  of  the  first  event, 
tne  same  distribution  F,  and  so  on.  So  wnen  a  event  occurs, 
we  say  that  a  renewal  has  taken  place. 

An  example  of  a  renewal  process,  let  us  suppose  that  we 
have  an  infinite  supply  of  lightbulbs  whose  lifetimes  are 
i.i.d,.  Suppose  also  that  we  use  a  single  lightbulb  at  a 
time  and  when  it  fails  we  immediately  replace  it  with  a  new 
one.  Under  these  conditions,  (N(t),  t>0  }  is  a  renewal 

process  when  N(t)  represents  the  number  of  lightbulbs  that 

have  failed  by  time  t.  Accordingly,  if  we  assume  that  this 
failure  is  a  kill  by  a  single  firer  -  X  or  Y  then  killing 
process  also  can  be  considered  as  a  renewal  process.  Thus, 
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renewal  (or  event  or  arrival)  will  be  synonymous  witn  <ill 

or  casualty  occurence  in  the  rest  of  this  thesis.  :’e  now 

assume  that  casualty  cccure nee  is  a  kind  of  renewal  process. 

If,  in  particular,  the  distribution  of  interarrival  time  of 

kill  is  exponential  with  probability  distribution  function 
x 

e  ,  then  this  renewal  process  turns  out  to  be  a  Poisson 
process  with  rate,  .  £  Bef  .  9] 

Thus  we  may  assume  tna t  X  casualty  occurence  due  to  a 
single  Y  firer  who  uses  the  same  acquisition  mode  to  engage 
each  X  target  can  be  explained  by  ordinary  renewal  process 
because  we  assume  that  times  between  casualties  are  all 
i.i.d.  However,  we  often  consider  a  counting  process  for 
which  the  first  interarrival  time  has  a  different  distribu¬ 
tion  than  the  remaining  ones.  So,  we  may  think  of  a  Y  firer 
who  at  first  uses  serial  acquisition  mode  to  engage  an  X 
target,  hut  once  he  has  finished  engaging  an  X  target,  then 
he  may  use  parallel  mode  to  engage  tne  rest  of  the  targets. 
So,  in  this  case,  we  can  assume  that  interarrival  time  of  X 
casualty  has  different  distribution  between  the  first  and 
the  rest. 

Formally,  let  {X  ,i>1)  be  a  sequence  of  independent  r.v. 

i 

with  X  having  distribution  G,  and  (X  ,i>2)  having  distritu- 
1  i 

tion  F  which  is  different  from  G .  Then  counting  process 

( t) ,  t2)  is  said  to  be  a  delayed  renewal  process.  So,  we 
can  think  of  using  delayed  renewal  process  for  different 
acquisition  mode,  but  here  in  this  thesis,  we  will  only  deal 
with  the  case  of  ordinary  renewal  process  for  simplicity. 

Xe  now  define  for  the  renewal  process  a  number  of  asso¬ 
ciated  random  variables  whose  study  is  the  objective  of  this 
thesis.  Table  IV  and  Figure  3. 1  give  the  intuitive  interpre¬ 
tation  of  the  random  variables  obtained  in  renewal  process. 

Here  we  can  interpret  each  random  variable  as  following: 

•  X  is  the  time  between  the  (n-l)st  and  n-th  X  force  kill 
n 
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TABLE  IV 

Random  Variables  in  Renewal  Process 


r.  v  interpretation 


X  time  between  the  (n- 1)  st  and  n-th  renewal 

n 

i-e.  n-th  interarrival  time 

S  time  ap  to  n-th  renewal 

n 

N  (t )  total  number  of  renewals  in  (0,  t] 

M  (t)  expected  number  of  renewals  in  (0,  t] 

known  as  renewal  function 


Z(t)  length  of  time  measured  from  the  last 

renewal  until  a  given  time  t  ,  called  as 
age  or  backward  recurrence  time 

Y (t)  length  of  time  measured  from  time 

t  until  the  next  renewal,  called  as 
excess  life  or  forward  recurrence  time 

S(t)  length  of  renewal  time  at  time  t, 

called  as  spread 


by  a  single  Y  firer. 

•  S  is  the  time  up  to  n-th  kill. 

n 

•  N  ( t)  is  the  total  number  of  kills  in  (0,t]. 
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S  (t) 


> 


<  — Z(t)  --> 


< — Y  (t)  — > 


x 


time 
-  x — 


s  s  s  t  s 

1  2  N(t)  N  (t)  +1 


Figure  3.1  Benewal  Process. 

•  M  (t)  is  average  number  oi  kills  in  (0,t]. 

•  Z  (t)  is  the  time  between  the  last  kill  and  a  given  time  t. 

•  Y  (t)  is  the  length  of  time  from  time  t  up  to  the 

next  kill. 

•  S  (t)  is  tae  length  of  lifetime  of  X  force  soldier  at 

time  t. 

In  this  thesis,  we  will  only  deal  with  X  force  casualty 
by  Y  firer  to  avoid  the  confusion  in  notation. 

C.  DISTEIB 0TIO N  OF  TOTAL  H0332B  OF  CASUALTIES 

Suppose  that  for  an  ordinary  renewal  process  a  single  Y 
firer  uses  the  parallel  acquisition  mode  to  engage  passive  X 
targets  from  tae  beginning  of  combat,  when  we  assume  that 
targets  are  already  acquired.  Then,  the  distribution  of 
N  (t) ,  the  total  number  of  X  force  casualties  by  a  single  Y 
firer,  can  be  obtained  at  least  in  theory  by  first  noting 
the  relationship  that  the  number  of  kills  by  time  t  is 
greater  than  or  equal  to  n  if  and  only  if  the  n-th  kill 
occurs  before  or  at  time  t. 

Formally 


/ 
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(3.1) 


N (t) >n  <==>  S  <t 
n 

From  3.1/  we  obtain 

?{N(t)=n}  =  P  £N ( t)  >n}  -  ?{N(t)>a+1) 

*  ?  £S  <t}  -  ?£S  <t} 

n  n+1 

Now,  since  the  random  variable  X  *s,(i£l),  are  independent 

i  n 

and  have  a  common  distribution  F,  it  follows  that  S  =  X 

n  i=1  i 

is  distributed  as  F  ,  which  is  n-fold  convolution  of  F  with 

□ 

itself  [Hef.  10].  Here,  F  is  tae  prooaoility  distribution 
of  time  between  casualties.  Therefore  we  obtain 

?£N(t)=n]  =  F  (t)  -  F  (t)  witn  ?  (t)  =  1 
n  n+1  0 

and 

P  £N  (t)  <n}  =  ?  £S  >t}  =  1  -  F  (t) 

n  n  (3.2) 

Therefore  the  probability  distribution  of  (t)  ,  which  is 

total  number  of  casualties  by  single  Y  rirer,  can  be 

obtained  explicitly  fcr  all  n. 

The  simplest  special  case  of  equation  3.2  is  obtained  by 

taking  the  renewal  process  to  be  a  Poisson  process. 

Then  S  has  the  special  Erlang  distribution  with  n  stages 
n 

because  each  F  is  exponentially  distributed.  Also,  it  is 
known  from  the  previous  argument  tnat  S(t)  has  a  Poisson 
distribution  with  mean  t. 

D.  EXPECTED  BOBBER  OF  CASUALTIES 
1 .  Renewal  Func tion 

Expected  number  of  casualties,  E[B(t)],  which  is 
defined  to  be  the  mean  number  of  casualties  in  the  interval 
(Q,t],  is  known  as  the  mean  value  0£  renewal  function.  B(t). 
From  equation  3.  1,  the  expected  number  of  casualties  in  time 
t  can  be  obtained  as  following: 
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.1(t)  =  E£N(t)  ]  -  I  k  P£»(t)»k} 

K-  1 

N)  T>J 

=F  ^  P  {H  (t)  =k} 
c= 1 k=n 

=  Z.PCS  (t)>n) 
n=  1 

c« 

=  r  ?£5  <tj 
n=  1  n 

So 

=  ^L  F  (t)  (3.3) 

e=  1  r. 


But  hare  it  is  generally  very  difficult  to  find  H  (t) 
directly  from  the  distribution  F.  Thus,  we  will  employ  the 
Laplace  Stieltzes  Transform  {L .  S.  T, )  [Ref.  9]  to  find  more 
conveniently  an  expression  for  M(t).  He  generally  need 
three  steps  to  find  M  (t)  . 


Step  1: 

From  the  definition  of  L.S.T.,  we  know  taat  L.S.T. 

of  a  (t)  is 

*  1  -st 

M  (s)  =  q  e  da  (t) 

and  also 

F(s)  e  dF(t) 


Step  2: 

Then  from  equation  3.3 

*  a-  * 

M(S)  ®  —  F  (S) 

n=  1  n 
$2  *  n 

*  —  *  F  (s)  ]  for  s>0 

n=  1 

The  above  is  a  sum  of  geometric  series  and  can  be  written  as 
follows: 


:i*s>  * 


As) 


1-F(s) 


(3.  4) 
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or  equivalently 


As) 


* 

M(S) 


* 

1+M(S) 


Step  3: 

Since  the  L.S.X.  of  a  function  uniquely  determines 
the  function,  the  renewal  function  determines  the  distribu¬ 
tion  of  the  kill  time  and  therefore  determines  the  prob¬ 
ability  law  of  the  casualty  occurence.  Specifically  there  is 
a  one-to-one  correspondence  between  interarrival  time 
distribution  F,  and  renewal  function,  d  (t)  .  So  if  we  have 
the  distribution  ?,  then  we  can  find  d(s)  from  tne  equation 
3.4.  But  to  find  3(t),we  have  to  transform  a  (s)  to  M  (t)  by 
inverse  application  of  L.S.  T. 

let  us  try  the  example.  Remember  that  this  example 
will  be  applied  all  the  way  through  this  chapter. 


Example  3. 1 

Suppose  that  each  inter  arrival  time  to  hill  X .  ,  (i>  1 )  has 
exponential  distribution  witn  .  Inen  L.S.T.  of  F  is 

*  -st 

?(s)  =  e  dF(t) 


/ 

f 

-  0 

'SO 


-st  .  -' t 

e  /•.  e  dt 


-  t 


dt 


\  +s 


-  (Vs)  t 

(Vs)  e  dt 


TVs 


if  s>- TV 


Then  from  equation  3.4 
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# 


Ms) 


* 

FjS) 


1-F(s) 


JV£.fL _ 


where  1/s  was  an  L.S.T.  of  F(t),  which  had  value  of  t.3  So 
M  (t)  =  /  t  ♦  c  (some  constant)  ,  but  we  know  that  d  (0)  =  0, 
because  there  would  be  no  casualty  at  time  0.  So,  finally 
M(t)  =  At. 

Let  us  be  more  specific  about  this  example.  Suppose 
that  interarrival  times  of  casualties  are  i.  i.  d.  and  expo¬ 
nentially  distributed  with  mean  30  minutes.  Then 

a)  What  is  the  distribution  of  total  X  casualties  at  5 
hours  later  ? 

b)  What  is  the  probability  that  3rd  X  kill  was  occurei  at 
time  t  =  4  ? 

c)  What  is  the  probability  that  at  time  t  =  2,  no  X  casu¬ 
alty  occured  ? 

To  answer  these  questions  we  know  tnat  3 (t)  *  2t  from  above 
example  because  =  2,  then  nuaner  of  X  casualties  in  time  t 
has  Poisscn  distribution  with  mean  2t. 


3Let  F(t)  =  t  for  t>0.  Ihen  L.S 

*  ( 50  -st  ^  -st 

F  (t)  -  \  e  d?  (t)  *\  e  Idt 

=  1/s  for  s>0 


of  F(t)  is 

1  -st 
- e 

s  •'  0 
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Then 


-10  n 
e  ( 10) 

a)  P{N  (5)  =n }  * - -  — 

n  ! 


3  -8 

8  e 

t)  P {N  (4)  =3}  - - =  .  0236  and 


0  -4 

4  e 

c)  P  {N  (2)  =0}  * - -  =  .0183 


2.  Benewal  Enua tion 

Now  aa  integral  equation  for  d(t)  may  he  obtained  by 
conditioning  on  the  time  of  the  first  casualty  occurence. 
Op on  doing  so,  we  obtain 

(70 

H  (t)  =  2:  N  ( t)  jX^  x]  if  (x)  (3.5) 

however 


E[N(t)|Xi*  x]  = 


•;  0 

V  1  ♦  H ( t-x) 


if  x>t 
if  x<t 


(3.6) 


for  if  the  first  casualty  occurs  at  time  x,  x<t,  then  from 
this  point  cn  the  process  starts  over  again,  and  thus  the 
expected  number  of  casualties  in  (o,t]  is  just  1  plus  the 
expected  number  of  casualties  in  time  t-x  from  the  beginning 
of  an  equivalent  renewal  process.  Putting  equation  3.t>  in 
3.5  yields 

•t 

a(t)  =  \  :  1*3(t-X)  ]dF  (x  ) 

=  F(t)  ♦  .  J  H  (t-x)  dF  (x)  (3.7) 
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where  equation  3.7  is  known  as  the  renewal  equation  and  r  is 
a  known  function  and  M  is  an  unknown  function  to  be  deter¬ 
mined  as  a  solution  to  the  integral  equation  3.7.  F.oss 
provided  the  solution  for  equation  3.7  [Bef.  11]  as 
following  and  the  solution  is  accepted  true  for  every  case. 

Expected  nusber  of  casualties  up  to  tiae  t  is 

a(t)  *  F(t)  ♦  {*  P(t-x)  dH(x)  (3.8) 

t  » 

£>■5 

where  H  (x)  =  H  F  (x) 
n=1  n 

So  it  allows  to  find  X{t)  using  only  the  known  distribution 
function  F,  which  is  the  distribution  of  time  between  casualt¬ 
ies  ,  thus  avoiding  finding  individual  F  . 

n 

3 .  Renewal  Density 

Now  let's  define  the  renewal  density  as  m(t)  = 
dd(t)/dt.  Remember  ttfat  ii(t)  is  the  expected  number  of  casu¬ 
alties  up  to  time  t.  So  a(t)  may  be  regarded  as  the  expected 
number  of  casualties  per  unit  tiae;  more  precisely,  the 
average  number  of  casualties  during  the  time  interval  from  t 
to  t+h  is 

/t+h 

M(t*h)  -  H(t)  =  m  (x)  d x  (3.9) 

Differentiating  the  renewal  equation,  we  obtain  the  renewal 
equation  for  the  density  m(t) 

a  (t)  *  f  (t)  ♦  \  f  (t-x)  m  (x)  dx  for  t>0  (3.  10) 

So, this  knowledge  of  the  density  f  of  interarrival  time  of 
casualties  is  sufficient  to  determine  the  density  m(x)  of 
casualty  occurence. 
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E.  LIMIT  BEHiVIOOfi  OF  CASOALTX  OCCUEESCE 

Be  know  that  total  number  of  casualties,  N  (t)  ,  go  to 
infinite  as  tine  goes  to  infinite,  when  we  assume  that  force 
size  is  infinite.  But  it  would  be  useful  to  know  at  what 
rate  H(t)  approaches  infinity.  That  is,  we  would  like  to  be 
able  to  say  something  about  limit  of  N(t)/t  as  time  goes  to 
infinite. 

From  the  definition  of  r.v.  we  know  that 

s  <  t  <  s 

N  (t)  N  (t)  ♦  1 

wnere  S  must  be  the  time  of  the  casualty  after  time  t. 

N  (t)  *1 

By  refer  to  [Bef.  10]  we  know  that 
N  (t)  1 

-  — >  -  as  t  — >  (3.  11) 

t  u 

/ 

where  u  =  E[ X  ]  =  [1  -  F(x)  ]dx  ,  which  is  the 

i  0 

expected  time  between  casualties. 

Here  the  function  1/a  is  often  called  as  the  rate  of  the 
renewal  process.  Thus  the  average  number  of  casualties  per 
unit  time  converges  to  1/u  .  Then,  how  aoout  the  expected 
average  number  of  casualties  per  unit  time?  Is  it  true  that 
K(t)/t  also  converges  to  1/u  as  t  -->  ?  This  result,  known 

as  the  elementary  renewal  theorem.  will  be  stated  without 
proof  [Hef.  11]. 

Elementary  B enewal  Theorem 

H  (t)  1 

-  - > - as  t  -->  Cv  (3.12) 

t  u 

Thus  whatever  the  distribution  F  is,  the  average  number  of 
casaalties  o£  tg  time  t  is  approximately  (for  large  t) , 
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,1 


t 

a 


(3. 13) 


H(t) 


Aiso,  the  average  number  of  casualties  in  the  interval 
Ct,  ilh)  for  h>0  is  approximated  ror  large  t  as  following. 


h  . 

8  (t+h)  -  fl{t)  - -  as  t  — >  CO 

a 

for  any  distribution  F  with  mean  j. 


(3.  14) 


Here  eguation  (3.14)  is  known  as  alacxwell 's  theorem  when  F 
is  not  lattice. 

From  the  example  3.1,  where  casualty  occurence  time  is 
exponentially  distributed  with  mean  30  minutes,  the  average 
number  of  X  casualties  up  to  time  t  is  approximated  by  t/u, 
which  is  2t.  Then  after  2  hours  of  combat,  we  can  assume 
that  the  average  number  of  X  casaalties  by  single  X  firer 
will  be  2t,  which  is  4  casualties. 

Another  limiting  result  which  may  prove  to  be  useful  in 
renewal  process  also  concerns  the  r.  v.  N  (t)  .  It  is  shown 
[Ref.  12]  that  N  (t)  has  an  asymptotically  normal  distribu¬ 
tion  with  mean  t/u  and  variance  t  T~£ u5 .  Thus 


—  .  v  -  t/u  N 

lim  P[M(t)<y}  - .aia-----  (3.15) 

t->Cu  —  •  j  _2  3 

\\  <T  t/u 

Suppose  of  the  same  example  3.1  that  casualty  occurence  mean 
time  is  30  minutes.  Then  u  =  0. 5  =  S[  T  ]  =  1/A  .  So  A  =  2. 
Also  Var[TJ  *  1/x1  -  1/4  =  0.25  hr"  ,  which  is  900  min1.  here. 
Let's  assume  that  we  had  a  10  hours  of  combat  action.  Then 
expected  number  of  X  casualties  is  £[N(1v))  ]  =  t/u  =  10/0.5  = 
20  and  var[N(10)  ]  =  t  jVus  =  10x0.  25/ (0. 5)3=  20. 

Suppose  of  the  guestion  that  how  many  X  soldiers  will  be 
required  in  order  to  with  probability  of  95%  that  at  least 
more  than  one  X  soldiers  will  survive  after  10  hours  of 
combat  ?  This  question  may  be  answered  as  following. 
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0.95  <  P  {H (10)  <x} 


x  -^20  >  -p  '  z  -  20  ' 
— \  \~TT  —  a. 4 72 


Now  from  the  normal  table,  0.95  percentile  of  normal  distri¬ 
bution  is  1.645.  So  1.645< (x  -  20) /4. 472,  where  the  answer 
turns  out  to  be  *>27. 36  .  So  we  can  say  that  X  force  needs 
at  least  28  soldiers  as  an  initial  force  size  to  te  9SZ  sure 
that  at  least  more  than  one  X  soldier  can  survive  after  10 
hours  cf  comba  t. 

There  also  exist  a  alternative  normal  approximati on. 
when  we  assume  that  for  integer  n, 

?{N(t)<n}  =  P  (S  >t)  where  S  =  x 

n  n  i=1  i 

from  equation  3. 1. 

Then  P{S  >t}  =  1  -  P  {S  <t}  =  1  -  3  /-- — ----  (3.  1o) 

n  n  <\ — i  n 

So  now  we  assume  that  S  has  normal  distribution  with  mean 

n 

nu  and  variance  n-f  .  If  we  apply  equation  3.16  to  previous 
example, 

i  10  -  0.  5x  \ 

0.  95<?{N  (10)  <x]  *  1  -  £.'■ - r - 

>•  0.  5  -i  x 


-  ,  10  -  0.5x  .  10  -  0.5x 

Then  ^  : - <  0.0  5  or  - - -  < 

~  0.5  iPx  i  0.5.)  x 


1. 645 


If  we  solve  for  x,  it  turns  out  to  be  x>28.83.  So  answer  is 
almost  same  when  we  assume  some  rounding  errors. 


F.  AGE  AND  EXCESS  LIFE  DIS TRIBUTION  OF  CASUALTY  OCCUHENCE 

As  we  mentioned  earlier  in  this  chapter,  we  can  also 
consider  such  random  variables  as  Y(t),  which  is  the  time 
from  certain  time  t  until  next  casualty  occurence. 
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vaich  means  residual  lifetime  of  a  soldier  in  combat  from 
certain  time  t  if  he  is  still  alive,  ana  2  (t)  ,  vhici :  is  ca« 
time  from  t  since  the  last  casualty  occurence  (see  Figure 
3.2). 


s 

H  (t) 

S 

N(t)*1 

*  •****-*  *X“  1  "  ^  ■ 

1 

2(t) 

Y(t) 

S(t) 

Figure  3.2  Age  and  Excess  Life  of  Casualty  Occurence. 


This  means 
Y(t)  -  S 


S  (t )  ♦  1 

Z  (t)  =  t  -  s 


-  t 


S  (t) 


where  Y  (t)  is  called  as  the  excess  life  or  forward  recur¬ 
rence  time  and  Z  (t)  is  called  as  ace  or  pack  ward  recurrence 
time.  Also 


S(t)  =  Z(t)  *  Y(t) 


=  s  -  t  ♦  t  -  s 

N(t)*1  N  (t) 


s  s  —  s 

N  (t)  ♦  1  N(t) 


i=1 


X  -  ¥» 


X 

i=1  i 


X 


.V(t>  +1 


(3.  17) 
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Here  S(t)  is  often  called  as  spread  and  S(t)  is  also  r.v. 

because  X  is  r.v.  too. 

N  { t)  1 

Suppose  of  the  Poisson  process  that  X.  is  exponential, 
then  I[S(t)  ]  =  E[¥  (t)  ]  ♦  E[2(t)]  =  1/*  ♦  1/,\X*  2/N  =  22£ X  1. 
So,  S  (t)  can  be  assumed  as  the  length  of  the  lifetime  of  tne 
1  force  soldier  which  is  in  combat  at  time  t. 


IV.  THE  SOPEBPOSITIOH  OF  CASH ALT I  OCCDBENCES 

A.  IHTSCDOCTIOH 

In  the  previous  chapter,  ve  have  considered  in  detail 
the  theory  of  casualty  occurence  through  the  ordinary 
renewal  process  of  a  single  Y  force  soldier.  We  deal  in  this 
chapter  more  briefly  with  the  superposition  of  several 
renewal  processes,  which  is  total  X  force  casualty  occurence 
behaviour.  When  we  assume  multiple  Y  force  soldiers  who 
shoot  at  X  force  soldiers  independently  of  each  other,  we 
may  assume  that  all  these  Y  soldiers  kill  X  soldiers  by 
independent  casualty  processes.  So,  we  can  say  that  there 
occur  multiple  independent  ordinary  renewal  processes, 
simultaneous  in  time. 


•x  — 


process  1  | 


process 

process 

pooled 

output 


Figure  4.1  Superposition  of  Casualty  Occurences. 
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Suppose  that  we  have  '  p'  independent  casualty  processes 

in  operation  simultaneously ,  ail  wita  tne  same  P.D.F.  of 

casualty  occurence  time.  Consider  tne  sequence  of  casualty 

occurences  formed  by  pooling  the  individual  processes.  Fig 

4. 1  illustrates  the  special  case,  when  p  =  3.  Formally,  we 

can  say  that  if  {N  (t),  i>0 } ,  for  i  =  1,2,  ...  ,  k,  are 

mde- 

i  k  pendent  renewal 

processes  and  N  (t)  =  a  (t)  ,  then  (N  (t),taO} 

p  i=1  i  p  is 

called  pooled  process  which  is  total  X  force  casualty  occu- 
-  rence  behaviour  by  k  7  force 

soldiers. 


B.  SCUE  GENERAL  PSOPEBTIE S 

1 .  Normal  Approximation  of  Pooled  Output 

From  casualty  occurence  processes  in  the  previous 
chapter,  we  assumed  that  there  exist  a  mean  casualty  occure- 

2 

nee  time.  Sc  let  u  and  r  be  the  mean  casualty  occurence 

i  i 

time  and  corresponding  variance  by  tae  i-th  Y  force  soldier 
and  let  u(p)  and  ^(p)  be  the  mean  casualty  occurence  time 
and  corresponding  variance  of  the  total  casualty  occurence  by 
k  of  Y  force  soldiers.  Then  the  question  is  how  are  these  two 
values  rela ted  each  other  ? 

From  previous  study,  we  know  that  1)  N'  (t)  is  approxima- 

i 

2  3 

tely  normal  with  mean  t/u  and  variance,  n  t/u  .  Therefore 

i  '  i  i 

we  can  say  that  M  (t)  is  also  approximately  normal  witn  mean 
k  ?  k  2  3 

2 T  t/ui  and  variance  (J~  t/u  .  2)  Also  assume  now  that 

i=  1  i=1  i  i 

either  tne  pooled  process  is  a  renewal  process  or  at  least 
the  analogous  result  holds,  then  approximately  N  ( t)  is  ncra- 

2  3  P 

al  with  mean  t/u  (p)  and  variance  <j  (p)  t/  u  (p). 

So,  from  1)  and  2)  ,  we  must  have 


t  k  1 

-  =  t  r  —  - 

u  (p)  i=  1  u 

i 


(4.  1) 
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k  1  -1 

then  a  { p)  =  ( 1C - )  which  is  the  harmonic  mean  of  tne  u  , 

i=1  u  i 


2  2 

(p)  k  i 

and  - -  t  =  t  n  - 

3  i  =  1  3 

u  (P)  a. 

i 


k  Mi 


k  1  -3  k 


then  r~(F)  =  u  (p)  - -  =  (.21—  ) 


i  =  1  3 

a 

i 


i  =  1  u 


1=1  3 

a 

i 


(4.2; 


Also  if  we  assume  that  u  =  u  and 


Ji 


for  all  i  =  1,2,  ..,k 


which  means  all  casualty  occurences  are  independent  and  identi¬ 
cal,  then  mean  casualty  occurence  time  and  corresponding  vari¬ 
ance  of  the  total  casualty  behaviour  by  k  of  Y  force  are  as 
follows; 


'  1  -1  u 

u (F)  =  (  -  k  )  =  - 

u  k 


2  _  2 

2  1  -3  T  T  J~  2 

T  (?)  =  (  -  k  )  —  k  =---=(-  ; 

u  3  2  n 

u  k 


(4.  3) 


Also  total  number  of  expected  casualties,  E[  N  (t)  ],  can 
be  obtained  here  by  pooling  individual  casualty  occurences 
as  following; 

EC*  it)  ]  *  £  E£».  (tj  ]  =  k  E[H.  (t)J  =  k  M  (t)  (4.4) 

p  1=1  i  i 

where  a  (t)  is  expected  number  of  casualty  from 
single  fixer,  which  is  t/u  for  large  t. 
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let's  look  at  an  example  that  deals  with  pooling  of 
individual  renewal  processes. 

Example  4.  1 

Suppose  of  X  force  artillery  shells  from  several 
sources  are  bombording  the  same  X  target.  Each  source  hurls 
shells  at  the  target  at  a  rate  of  40  saells  per  hour.  Assume 
that  the  intar  arrival  tiaes  of  shells  from  each  single 
source  are  uniformly  distributed  over  an  interval  (0,  a]. 
Then  the  question  is  1)  What  is  the  probability  that  less 
than  825  shells  will  be  nurled  at  the  target  from  a  single 
source  during  a  20  hour  period  ana  2)  If  there  are  5  sources 
from  which  such  attacks  are  launcned,  what  is  the  prob¬ 
ability  that  more  than  4125  (  which  is  5  times  825  )  shells 
are  nurled  at  the  target  in  a  20  hour  period  ? 

To  answer  the  first  question,  let  ^  =  40  shells/hour 
=  2/3  shells/min.  Then  u  =  E[  X  ]=  1/~  =  3/2  min,  when  we 
assume  pcisscn  process.  3ut  interarrivai  times  are  uniformly 
distributed  over  an  interval  (0,a].  So,  u  =  E[X]  =  (0  + 
aj/2=  a/2,  where  a  =  3  minutes.  Then  Var[X]  =  (a  -  0)‘/12  = 
3/4  min2"  .  Now  let  N  (t)  =  number  of  shells  hurled  in  (0,  t] 
from  a  single  source.  Since  K(t)  us  normal  with  mean  t/u  and 
variance  tr/u3  from  equation  3.15,  for  t  -  20  hours  =  1200 
minutes,  the  probability  turns  out  to  be  as  following; 


325  - 


?  {N  (1200  )  0  2  5}  =  5. 


—  '  825  -  800 
~  16.33 


“  .  120  0  -  3/4  —  16.33 

rivi” 

=  3(1.53)  =  0.9370 
and  for  second  question. 

Let  N.(t)  =  number  of  shells  hurled  in  (0  ,  t]  by  i-th  source 
where  i  =  1,2,  ...  ,  5  .  and 


5 

N  (t)  =  72  N  (t)  =  total  number  of  shells  hurled  in  10,  t 
p  i=1  i 


40 


by  5  sources. 

Then  {M  (t)  ,  t>0}  is  the  pooled  process  for  which  mean 
p 

lr.terarnva  1  time  is, 

u  3/2  3 

u  (p)  =  -  =  - -  —  minutes. 

k  5  10 

and  corresponding  variance  is 
(f  (?)  =  (  aiQ 

Since  N  (t)  is  normal  with  mean  t/u  (p)  and  variance  V  (F)  t/</  (p) 
for  t  =  *1200  minutes, 

1200 

4  125 - x 

:c  I . .s  'All..  '■ 

=-  1200  3/100  — '  36.51 

- -  e 

u  3 

\  (3/10) 

=  5(3.4  2)  =  .9  997 
So,  P  (1200)  >4  125}  =  0.0003 


?  (N  ( 1200)  <  4  1  25}  = 

P 


2.  Poisson  Process  of  Pooled  Output 

Now  let's  assume  that  individual  casualty  processes  are 

probabilist ically  identical  witn  common  underlying  C.D.F.  of 

casualty  occurence  time  r(x).  Let  F  (x)  be  the  C.D.F  of 

p 

casualty  occurence  time  of  pooled  process.  Then  the  question 

is  how  is  F  (x)  related  to  F  (x)  ?  This  question  can  be 

answered  byrusing  the  idea  of  excess  life.  Let  Y  (t)  be  the 

i 

excess  life  of  an  X  force  soldier  and  Y  (t)  be  excess  life 

of  total  X  force.  Then  we  know  that  Y  (?)  *  min[  Y  jt)  ,  Y  (t^ 

_  ,  Y  (t)  j  (see  Fig  5.  1)  T 

k 

Then 


4  1 
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Figure  4.2  Superposition  of  Excess  Life. 


k 

P(Y  (t)  >y}  =  ?{Y.  (t)  >y} 

P  i 


(4.  5) 


From  previous  study,  we  know  that  P.D.F  of  E.P.  for  individual 
renewal  processes  is  ( 1/u)  [  1  -  F  (x)  ]  and  then  we  can  say  that 

1 

—  [1  -  F  (x)  ]  is  P.D.F  of  Z.  D.  for  the  pooled  casualty 
u(p)  p 

process.  Thus, 


?{Y(t)>y} 

P{Y  (t)  >y} 
P 


,-CO 


1 


,  —  £  1  —  F  (x)  ]dx  then , 

y  u 


X) 


1 


*  \  ----  :  i  -  f  (x)  ]dx 

■y  «(p)  p 

1  50  k 

=[  -  \  [  1  -  FU)  Jdx] 

u  y 


from  equation 


4.  5 


If  we  differentiate  with  respect  to  y 


1  1  p  k- 1  1 

-  £  1  -  F  (y)  ]  «  k[  -  .  [  1  -  ?  (x)  ]dx  ]  (-  -  ) ;  1  -  r  (v) 

u(p)  p  ay  u 

Here  we  krow  that  u(p)  =  k/u,  so  a's  are  all  caret! id  cat. 
Then, 


C.D.F.  of  pooled  output  is 

i  k- 1 

F  (y)  -  1  -  £  -  \  [  1  -  F  (x)  ]dx  ]  [  1  -  F  (y)  ]  (4.6) 

P  u  y 


For  example,  think  of  the  exponential  distribution  of 

-  X  x 

casualty  occurence.  Then  F(x)  *  1  -  e  ,  and  u  =  “\/K  . 

Then  C.D.F  of  pooled  casualty  occurence  is 

->x  k-1  -\y 

F  (y)  =  1  -  £  \  [  e  ]d  x  ]  e 

p  y 


-Ay  k-1  ->y 
=  1  -  (e  )  e 


*  1  -  e 


-\ky 


for  v>0 


which  follows  that  pooled  output  is  a  Poisson  process  with 
rate  k  where  k  is  the  number  of  I  force  combatants. 

Here,  probably  the  most  interesting  properties  of 
the  pooled  output  refer  to  the  'local'  behaviour  where  k  is 
large.  Khintchine  [  Ref .  13]  has  proposed  that  in  the  limit 
the  numbers  of  renewals  in  non-overlapping  intervals  follow 
independent  Poisson  distributions,  thus  showing  that  in  the 
limit  the  cooled  output  is  §  Poisson  process.  Also  his 
proof  does  not  reguire  each  interarrival  time  to  be  identi¬ 
cally  distributed.  In  attrition  analysis  in  combat,  it 
appears  to  be  the  same  way,  which  means  that  if  a  large 
number  of  independent  (and  identical)  casualty  occurences 
are  pooled,  then  the  total  force  casualty  occurence  is 
approximately  a  Poisson  process  [Hef.  9]  [Ref.  13].  Thus, 
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we  can  say  that  pooled  casualty  occurence  by  total  Y  coafa- 
tan t s  c^n  be  approximated  b  y  Poisson  process  if  we  as suae  a 
lar je  cumber  of  Y  forces. 

Thus,  finally,  from  the  oasic  theorem  of  Poisson 
process,  we  know  that  if  (N  (t)  ,t>0}  is  Poisson  process  with 
rate  >0,  then  the  interarrival  tines  {X.,i>1}  are  i.i.d. 
exponential  random  variable  with  rate  A  .  Thus,  from 
previous  results,  total  X  casualty  occurences  is  a  Poisson 
process  with  rate  ,  where  k  is  the  number  of  Y  force 
combatants,  and  interarrival  times  of  X  casualties  can  be 
considered  to  be  exponentia llv  distributed  with  parameter  k 


V.  CCflPABISON  21  BEMEHAL  PBOCESS  WITH  COM POT EE  SIMULATION 

In  this  chapter  we  will  loo*  at  the  results  of  seme 
typical  combat  attrition  process  (  Y  firers  with  different 
number  of  passive  X  targets)  obtained  through  the  theoret¬ 
ical  renewal  process  model  which  is  developed  in  this  thesis 
and  computer  simulation  using  random  numoer  generator. 

As  we  mentioned  earlier,  total  number  of  expected  casu¬ 
alties  by  n  Y  firers  turned  out  to  be  n  M(t)  for  large  t 
from  the  renewal  theory,  where  we  assumed  that  all  casualty 
occurences  were  independent  and  identical.  But  it  is  clear 
that  these  expected  casualties  can  not  be  bigger  than  the 
number  of  X  targets.  Now  we  will  compare  this  result  with 
simulated  output  to  see  whether  we  can  find  any  interesting 
facts  such  as  how  many  X  targets  will  ne  most  appropriate  to 
apply  renewal  theory  to  combat  attrition  ?  or  what  force 
ratio  will  be  most  applicable  to  renewal  process  ?,  etc. 

In  order  to  facilitate  the  simulations,  the  assumptions 
were  made  that  n  firers  shoot  independently  at  m  targets, 
and  that  the  times  to  hill  the  targets  are  uniformly 
distributed  between  [0,1].  When  a  target  is  killed,  all 
firers  shooting  at  it  are  assumed  to  select  a  new  target  at 
random  from  those  still  surviving. 

The  procedure  of  the  computer  simulation  is  as  follows 
(see  Appendix  A)  ; 

1.  We  generate  the  initial  kill  time  for  each  firer  by 
using  pseudo  random  number  generator,  where  we  assume 
that  every  firer’s  kill  time  is  uniformly  distributed 
between  [  0, 1  ]. 

2.  Each  target  is  indexed  by  integer  number. 

3.  We  need  to  decide  which  target  does  the  firer  snoot 
at  ?  The  target  selection  rule  is  such  that  first  we 
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generate  a  random  uniform  number  between  [0,1],  then 
we  take  the  integer  value  of  (random  number  times  the 
number  of  targets  )  plus  one,  which  designates  the 
index  of  target  that  is  assigned  to  firer  1.  Same 
procedures  are  repeated  until  all  the  firers  are 
assigned  a  target,  with  some  targets  possibly  hcving 
more  than  one  firer.  Thus  at  first  tine  step  random 
number  for  selecting  the  target  is  generated  n 
times,  because  we  have  to  assign  a  target  to  each 
firer.  These  targets  are  killed  eventually  as  time 
proceeds. 

Now,  the  target  that  will  be  killed  first  is  the  target 
which  is  assigned  to  the  firer  who  has  the  minimum  time  to 
kill  (  *  Tain)  among  the  firers,  where  this  time  to  kill  was 
already  generated.  So  we  now  know  which  target  is  killed. 

After  the  first  target  is  killed,  the  number  of  targets 
killed  is  registered  as  one  and  the  first  time  step  ends. 
Now,  in  the  second  time  step,,  the  new  kill  times  between 
£0,1]  are  generated  again  for  all  firers  who  were  shooting 
at  the  first  target.  The  firers  who  were  still  in  the 
process  of  killing  a  target  at  tae  end  of  the  first  time 
step  continue  to  engage  the  same  target  in  the  second  time 
step,  where  time  to  kill  in  the  second  step  is  old  kill  time 
{  generated  value  at  first  time  step  )  less  Tmin. 

Cnee  the  second  target  is  killed,  the  second  time  step 
ends  and  the  number  of  targets  killed  is  now  registered  as 
two.  Now  we  proceed  to  the  third  time  step,  where  same 
target  selection  rule  and  firing  rule  are  applied  until  all 
the  targets  are  killed. 

After  all  the  targets  are  killed,  each  time  step  length 
(  which  is  the  same  as  the  time  to  kill  the  target  in  each 
time  step  )  is  stored,  and  this  is  one  replication  (see 
Appendix  B) .  To  obtain  the  expected  value,  40  replications 
of  the  simulation  were  conducted  and  6  data  sets 
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(11=3,  m=6)  ,  (d  =3,m=  1  0) ,  (n=  3,  m=14)  ,  and  (n=1Q,  a=o)  ,  ( n=  10 
m=10)  ,  (is  1 0,  111=14),  were  examined.  After  40  replications 
the  data  were  evaluated  to  find  the  expected  number  of  kills 
through  order  statistics.  This  expected  numcer  of  kills  was 
plotted  against  time.  The  results  for  n=3  and  n=10  are 
shown  in  Fig  5.1. 

The  theoretical  results  can  be  obtained  through  two 
methods,  which  are  discrete  approximation  and  simulation. 
In  this  thesis,  we  used  discrete  approximation  based  on  the 
renewal  function  modelled  at  equation  3.7.  It’s  formula  is 
as  follows; 


M  (t)  =  i'L  1  ♦  n  (t-j)  ]  P . 

3  =  1  3 

where  P  represents  the  probability  distribution 

j 

of  kill  time. 

Thus,  expected  number  of  kills  obtained  from  the  renewal 
process  model  using  discrete  approximation  (see  appendix  C 
for  program  and  data,  where  left  column  refers  to  the  kill 
time  and  right  column  refers  to  the  expected  number  of  kills) 
is  alsc  plotted  against  time  in  Fig  5.1. 

It  was  found  through  both  cases  that  the  expected  number 
of  kills  by  the  renewal  process  model  may  only  be  justified 
through  the  assumption  of  an  infinite  number  of  targets 
because  as  m  gets  bigger,  the  expected  kill  value  approaches 
the  theoretical  result.  Also  it  was  discovered  that  at  the 
beginning  of  combat,  the  attrition  process  almost  followed 
the  theoretical  renewal  process  model. 
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VI.  CONCIOSIONS 


A.  SOflHARI 

Since  Lanchester  proposed  his  model,  many  studies  have 
tried  to  model  the  ccmbat  attrition  process  mathematically. 
Recently,  the  main  focus  has  been  on  calculation  of  numer¬ 
ical  values  for  Lanchester  attrition-rate  coefficient.  This 
numerical  determination  stands  at  tne  heart  of  casualty 
assessment  in  such  models.  There  are  three  reasons  for  the 
importance  of  this  attrition-rate  coefficient  methodology. 
First,  Lanchester- type  models  are  used  in  various  U.S.Army 
and  D.O.E.  planning  activities  more  widely  than  ever  before. 
Second,  Lanchester  attrition-rate  coefficient  is  a  basic 
element  of  any  Lanchester-type  combat  model,  and  that 
attrition-rate  coefficient  reflects  tne  effective  applica¬ 
tion  of  firepower.  Finally,  significant  new  developments 
have  occured  in  methodology  for  developing  more  tactically 
realistic  Lanchester  attrition-rate  coefficient  and  these 
important  results  have  not  been  accessible  to  a  very  wide 
audience.  In  particular,  a  new  approach  for  developing  more 
realistic  and  mathematical  Lanchester  attrition-rate  coeffi¬ 
cient,  that  of  computing  combat  attrition  by  using  a  renewal 
process,  was  proposed  by  Iaylor  recently.  This  thesis  began 
with  the  assumption  that  the  combat  attrition  process  was  a 
kind  of  renewal  process,  and  that  times  between  casualties 
were  i.i.d.  Thus,  considering  the  problem  from  the  stand¬ 
point  cf  renewal  theory,  many  new  random  varianles  were 
included  and  various  distribution  functions  studied  in  order 
to  understand  better  this  type  of  casualty  occurence 
paradigm. 


This  thesis  only  dealt  with  the  case  of  hono  j  ir.ecus 
forces  which  were  using  the  same  acquisition  mode  to  attrit 
the  targets,  which  were  passive  X  forces.  la us  we  could  use 
an  ordinary  renewal  process  to  study  the  behaviour  of  casu¬ 
alty  occurences.  Me  emphasized  the  probabilistic  function 
ox  each  r.v.,  so  once  we  found  the  distribution  of  interar¬ 
rival  times  of  casualties,  casualty  analysis  could  be  done 
by  plugging  distribution  functions  into  the  equation  devel¬ 
oped.  Finding  the  distribution  of  interarrival  tines  of 
casualties  was  the  critical  factor  in  continuing  the  casu¬ 
alty  analysis. 

The  final  conclusion  is  that  total  casualty  occurence  in 
coab^t  can  be  assumed  to  be  a  Poisson  process.  and  thus 
tiaes  between  casualties  ajre  exponentially  distributed  if  we 
ass  me  4  large  number  of  coabatants. 

Me  have  dealt  here  only  with  a  passive  X  force  (  i.e.  one 
that  does  not  shoot  back  )  ,  perhaps  an  unrealistic  assump¬ 
tion.  But  still  this  study  permits  better  understanding  of 
combat  attrition,  through  the  use  of  probability  theory, 
than  was  previously  available. 

B.  SUGGESTIONS  FOR  FURTHER  STUDY 

In  this  thesis,  we  have  only  dealt  with  passive 
targets.  A  more  complex  formulation  may  be  required  for  real 
combat,  since  enemy  targets  shoot  back  and  cause  firer 
attrition. 

In  addition,  a  delayed  renewal  process  could  be  applied 
for  the  case  where  different  acquisition  modes  are  used  for 
the  first  shooting  and  for  each  subsequent  shootings.  A 
reinforcement  of  the  forces  could  also  be  included. 
Casualty  analysis  could  then  be  done  using  the  theory  of 
alternating  renewal  processes,  once  we  know  the  rate  of 
reinforcement. 


Ag  gENDIX  A 

FORTRAN  PROGRAMMING  FOE  SIMULATING  THE  KI1L  TIME  (N  =  3) 


100 

105 


110 

10 


130 


125 


70 


135 

200 

300 


GO  10  130 


REAL  U1 (3)  -0  2  (3) -0  3,  U 4, THIN. TIME  (3)  . INTIKE (40 ,  14) 
INTEGER  JXGl  ,  F INDEX .L  ,  M,  N,  KTGT,  T&T  (\  4)  ,  RIG T  {14} 
DOUBLE  PRECISION  DSEED 
CSZED  =  123457. DO  C 
DO  20  0  K  =  1,4  0 
M  =  14 
N  =  3 
I  =  1  C 

CALL  3GUBS (DSEED, N,  U1)  C 
DO  100  I  =  1.N 
TIME  (I)  =  01(1) 

CONTINUE  C 
DO  105  J  =  1,M 
HTGI(J)  =  5 
CONTINUE  C 

DSSZD  =  DSEEC  +  23.  DO  C 
CALL  3GUBS  ( DSEED/N,  U2)  C 
DO  110  I  *1,1 

TGT(I)  =  I  NT  (U2  (I )  *  FLOAT  ( M)  )  «•  1 
CONTINUE  C 
CONTINUE 
TMIN  =  2. 

CO  130  I  *  1,N 

IF  (TMIN.  LE.  TIME  (I)  ) 

I  MIN  a  TIME  (I) 

KTGT  *  TGT  (I) 

FINOEX  *  I 
CONTINUE  C 

INTIME (K  ,L)  =  TMIN  C 

DO  125  J  *  1  ,M 

IFjRTGTjJ)  .  NE.  KTGT)  GO  TO  125 

RTGT  ( JTGT)  =  RTGT  (M) 

CONTINUE  C 
M  =  S  -  1 
1  =  L  ♦  1C 
DO  135  I  =  1,N 

IF(KTGr.  EQ.TGT(I)  )  GO  TO  70 
TIME  (I)  =  TIME  (I)  -  IM IN 

GO  TO  135 

DSEED  =  DSEED  +  1  2. DO 
CALL  GGUBS  (ESESD,  1,U3) 

TIME (I)  =  U3 
DSEED  =  DSZED  +  3 .  DO 
CALL  GGUBSlDSEED,  1.U4) 

T  GT  ( I  f  =  RTGI(INT  (J4*rL0AT(M)  )  *1) 

CONTINUE  C 
IF  (M.3E.  1)  GC  TO  10 
CONTINUE 

WHITE  (6,30  0)  (  (INTIME  (X,L)  ,  L=  1 ,  1  4)  ,  K=  1 , 4  0) 

STOP 

F0EMAT(1X,  14F5.2) 

END 


51 


APPENDIX  3 

SIMULATED  DATA  FOB  KILL  TIME  (N  =  3,  M  =  6) 


0.  260  7 

0.  5056 

0.2000 

0. 4392 

0.1057 

0.0649 

0.  3914 

0.  3860 

0.  16  59 

0.  3315 

0. 0303 

0.1646 

0.  1156 

0.  2987 

0.01 87 

0.391  1 

0.  0438 

0.4744 

0. 4545 

0.  0154 

0. 13  12 

0.2905 

0.1360 

0.3672 

0. 4133 

0.  1269 

0.0678 

0. 13o5 

0. 1 724 

0.0101 

0-  0834 

0.  0730 

0.2945 

0.0234 

0. 7305 

0.0620 

0. 4117 

0.  3926 

0.03  63 

0.0c97 

0.0277 

0.0151 

0. 0499 

0.  24  10 

0.3564 

0.0033 

0.  1442 

0.2979 

0. 463  7 

0.  2188 

0.  1503 

0.054b 

0.  5555 

0.0230 

0.7428 

C.  0736 

0.36  29 

0. 0375 

0.0319 

0.0322 

0.  232  6 

0.  3066 

0.2037 

0.2382 

0.  1050 

0 . 1466 

0.  1105 

0.  50  15 

0.23  16 

0.  1160 

0.2916 

0.4032 

0. 0943 

0.  1942 

0.2346 

0.0342 

0.2771 

0.0218 

0.  5372 

0.  01  04 

0.36  58 

0. 148  1 

0. 0724 

0.1435 

0. 0492 

0-  0554 

0.  1336 

0.  1750 

0.5712 

0.0609 

0. 2496 

C.0  2  27 

0. 11 02 

0.2093 

0.2394 

0.1294 

0.  453  0 

0.  0027 

0.1524 

0. 0436 

0.0132 

0 . 62  76 

0. 2916 

0.  2342 

0.0265 

0.5775 

0.301b 

0.1675 

0.  2299 

0.0756 

0.5040 

0.  1737 

0. 0637 

0.5797 

0.7527 

0.  0343 

0.1628 

0.0063 

0 . 3606 

0.0877 

0.  2201 

0.3  036 

0.1323 

0. 557  0 

0.0256 

0.3332 

0.609  7 

0.  0538 

0.05  16 

0.  333  1 

0.  1933 

0.0139 

0.7823 

0.  0719 

0.1470 

0.015b 

0. 074 3 

0.0118 

0.  1133 

0.  6589 

0. 0 1 44 

0.0712 

0.1216 

0.2234 

0. 0035 

0.  0602 

0.44  10 

0.  2291 

0.  1626 

0.2609 

0.  0174 

0.  56  6  9 

0. 1626 

0 . 20o  2 

0.1578 

0.0334 

Q.6032 

0.  00C3 

0.26  92 

0. 0024 

0.  0723 

0.2653 

0.3186 

0.  04  92 

0.3918 

0. 1603 

0.1857 

0 . 03  o  9 

0. 5564 

0.0105 

0.0674 

0.2037 

0.  018  4 

0.1233 

0. 5706 

0.  3  844 

0.3504 

0 .2022 

0. 0417 

0.1048 

0.  2102 

0.  4294 

0.11 55 

0. 0393 

0.  1705 

0 . 4 1  o  5 

0. 019  1 

0.  2929 

0.4283 

0.  2152 

0.  1427 

0.4954 

0.  1570 

0.0601 

0.3111 

0.  4439 

0.  0499 

0.0377 

0.  1931 

0.  46  58 

0.0251 

0.  15o8 

0.0375 

0.2395 

0.  1372 

0.41  83 

0.1055 

0.2174 

0.0335 

0.2340 

0. 0523 

0.  2950 

0.4356 

0. 0230 

0. 1 054 

0 . 0  9  2o 

0. 3024 

0.  5260 

0.06  69 

0.0079 

0.  1682 

0.3743 

0.  0645 

0.  15  04 

0.0088 

0.  0741 

0.  5531 

0.0228 

0.0220 

0.  43  50 

0.2320 

0.0207 

0. 2b34 

0.2070 

0.  128  5 

0.  1379 

0.1063 

0.0&95 

0.  0050 

0.5933 

IP  PSNDIX  c 

DISCRETE  APPROXIMA TIGN  FO  BENESA1  FUNCTION 


1 

2 

3 

4 

5 

6 

7 

8 


$  JCB 


203 

103 


REAL  FO  ,  M 1  (401)  ,  M  (41)  ,11X2  (41) 

pc  (=}  1 7/ 1 0  0  . 

DO  100  I  =  1*400 

IF  (I. 31.  1  00)  30  TO  230 
HI  |l+U  =  HI  (I)  ♦  (1.  ♦  XI  (I)  ) 
GO  10  100 

>1)  =  Ml  (I)  *  (1.  +  Ml  (I)  ) 
♦  SI  (1-10  0)  )  *  20 

CONTINUE 


10 

DO  11 0  J 

=  M1 

=  0.1*  (J- 1) 

11 

TIME  (J) 

12 

13 

113  CONTINUE 

Ml  (1  +  10*(J-1)  ) 

,*)  TIME  (J)  ,  M  (J) 

14 

15 

STOP 

16 

END  0 

SENTRY 

0.  0  0  0000  0 

0. 3000000 

0.  10  03030 

0.  1046219 

0.  2000000 

0. 2201895 

0.  300303  1 

0.3478478 

0.400000  1 

0.  4853622 

0.  5  0  30  00  1 

0. 644o293 

0.  600030  1 

0. 3166943 

0.  7000002 

1. 0067590 

0.  8000002 

1.  21  67 3  40 

0.  9000002 

1. 4486160 

1.  0  0  00000 

1.  7047900 

1.  09  999  9  0 

1.6733930 

1.  1  999990 

2. 0537170 

1.  3  000000 

2.  2452650 

1.  3939990 

2.  4373700 

1.  5300000 

2.  6341  510 

1.  6  0  CO  0  0  0 

2.  8344  850 

1.  6999990 

3.  0369570 

1.  3000000 

3.  2393290 

1.  8999990 

3.440^610 

2.  0000000 

3. 6377740 

2.  100000  0 

3. 8320340 

2.  1999990 

4.  023  1400 

2.  3  000  00  0 

4. 2256030 

2.  4000300 

4.  4239440 

2.  5  0  00  00  C 

4. 6227140 

2.  6000000 

4. 3215160 

2.  6  999  990 

5.  020)490 

2.  8000000 

5. 2181320 

2.  9  0000  0  0 

5.  41  57760 

3.  0000000 

5. 6132290 

3.  1000000 

5.  81  09340 

3.  2000000 

o. 3089100 

3.  3000000 

o. 2070170 

3.  4000000 

o. 4051470 

3.  5000000 

6. 6032300 

3.  6000000 

o. 30 12370 

3.  7000000 

6. 9991710 

3.  8000000 

7.  19  73  oOO 

3.  9000000 

7.  3549500 

4.  0000000 

7.  5923790 

PC 
*  PO 
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